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Abstract 

We sum up the next-to-next-to-leading logarithmic virtual electroweak corrections 
to the high energy asymptotics of the neutral current four-fermion processes for light 
fermions to all orders in the coupling constants using the evolution equation approach. 
From this all order result we derive finite order expressions through next-to-next-to 
leading order for the total cross section and various asymmetries. We observe an 
amazing cancellation between the sizable leading, next-to-leading and next-to-next-to- 
leading logarithmic contributions at TeV energies. 



PACS numbers: 12.38.Bx, 12.38.Cy, 12.15.Lk 



1 Introduction 



Experimental and theoretical studies of electroweak interactions have traditionally explored 
the range from very low energies, e.g. through parity violation in atoms, up to energies 
comparable to the masses of the W- and Z-bosons, e.g. at LEP or the Tevatron. The advent 
of multi-TeV-colliders like the LHC or a future linear electron-positron collider during the 
present decade will give access to a completely new energy domain. 

Once the characteristic energies s are far larger than the masses of the W- and Z-bosons, 
Mw,z, multiple soft and collinear gauge boson emission is kinematically possible. Conversely, 
exclusive reactions like electron-positron (or quark-antiquark) annihilation into a pair of 
fermions or gauge bosons will receive large negative corrections from virtual gauge boson 
emission. These double logarithmic "Sudakov" corrections PjJ |2] proportional to powers of 
g 2 In {s/My/ z ) are dominant at high energies and thus have to be controlled in higher orders 
to arrive at reliable predictions. 

The importance of large logarithmic corrections for electroweak reactions at high energies 
in one-loop approximation which may well amount to ten or even twenty percent was noticed 
already several years ago [3J Hj . The need for a resummation of higher orders of these double- 
logarithmic terms in the context of electroweak interactions was first emphasized in [3] which 
also contains a first discussion of this resummation. In particular, it was shown that the 
double logarithms do not depend on the details of the mechanism of the gauge boson mass 
generation. The issue is complicated by the appearance of massive (W, Z) and massless 
(7) gauge bosons in the SUl(2) x U(l) theory, which necessarily have to be treated on a 
different footing. A complete analysis of this problem in the double or leading logarithmic 
(LL) approximation by the systematic separation of soft (cj 7 < M) and hard (u y > M) 
photons was given in In two loops the results of this approach essentially based on the 
concept of infrared evolution equations have been confirmed by explicit calculations in [7j. 

The large coefficient in front of the single logarithmic term in the one-loop corrections to 
the electroweak amplitudes (see, e.g. jH]) suggests that sub leading terms play an important 
role also in higher orders, as long as realistic energies of order TeV are under consideration. 
Motivated by this observation a systematic evaluation of the next-to-leading logarithmic 
(NLL) terms for the neutral current massless four fermion process was performed in ref. 
Indeed one finds sizeable two-loop effects both for the total cross section, for the left-right and 
for the forward-backward asymmetry. A subclass of NLL corrections for general electroweak 
processes was subsequently evaluated in without, however, the very important angular 
dependent contributions. 

Various authors have also extracted the double and single logarithmic corrections from the 
complete one-loop calculations JUE]. The analysis for the general electroweak processes 
given in jT2] is in full agreement with [Hj, whereas a different prescription to separate the QED 
contribution is adopted in |TT]. Higher order heavy fermion mass effects on the asymptotic 
high energy behavior of the electroweak amplitudes were discussed in [T3|. The incomplete 
cancellation of the real and virtual electroweak double logarithmic corrections in the inclusive 
cross sections was investigated in [Ti] . 
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Following the approach of jH] which in turn is based on the investigations in the context of 
QCD [Ta[inilI71[IHllinil2niinil2ai231, this paper is devoted to the derivation of the next-to- 
next-to-leading logarithmic (NNLL) terms for the massless neutral current four-fermion cross 
sections. In Section 2 we present as a first step the NNLL form factor which describes the 
scattering amplitude in an external Abelian field for the SU(N) gauge theory. The derivation 
is based on the evolution equation derived in ref. [TBI EH [2D] • In Section 3 we then generalize 
the result to the four fermion process in the SU(N) gauge theory. After factoring off the 
collinear logarithms we use an evolution equation for the remaining amplitudes which is 
governed by an angular dependent soft anomalous dimension matrix [201 EU E3] • 

Finally we apply this result to electroweak processes in Section 4. To identify the pure 
QED infrared logarithms which are compensated by soft real photon radiation we combine 
the hard evolution equation which governs the dependence of the amplitudes on s with the 
infrared evolution equation j^J. The latter describes the dependence of the amplitude on a 
fictitious photon mass which serves as an infrared regulator and drops out after including the 
effect of the soft photon emission. The hard and infrared evolutions are matched by fixing 
the initial conditions at the scale M^z- Section 4 also contains a discussion of the numerical 
implications of our result. A brief summary and conclusions are given in Section 5. 



2 The Abelian form factor in the Sudakov limit 

Let us first consider the vector form factor which determines the fermion scattering amplitude 
in an external Abelian field for the SU(N) gauge model. In the Born approximation, 



•Fb = ipfahvrtipi) , (1) 

where p\ denotes the incoming and P2 the outgoing momentum. 

There are two "standard" regimes of the Sudakov limit s = {pi — P2) 2 — > —00: (i) on-shell 
massless fermions, p\ = p\ = and gauge bosons with a small non-zero mass M 2 <Q 2 [2] , or 
(ii) slightly off-shell fermions p\= p\ = —M 2 , and massless gauge boson pQ. Let us consider 
the first case and choose, for convenience, p\^ = (Q/2, 0, 0, =F<3/2) so that 2pip 2 = Q 2 = s. 
The asymptotic Q-dependence of the form factor in this limit is governed by the evolution 
equation [TH[ ITHj 







<91nQ : 



Q 2 dx 



M 2 X 



7 («(x)) + C(«(Q 2 ))+e(«(M 2 )) 



Its solution is 



f fQ 2 dr 

F = F (a(M 2 ))exp I — 
[Jap x 



x A T ' 

M2 — 7(«(x')) + C(«(a:))+^(M 2 )) 



(2) 



•F b ■ 



(3) 



The LL approximation includes all the terms of the form a n log 2n (Q 2 /M 2 ) and is determined 
by the one-loop value of 7(a). The NLL approximation includes all the terms of the form 
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a n \og 2n ~ m (Q 2 / M 2 ) with m = 0, 1. This requires the one-loop values of 7(a) and £(a) + £(a) 
and using the one-loop running of a in 7(a). The NNLL approximation includes all the terms 
of the form a n log 2n ~ m (Q 2 / M 2 ) with m — 0, 1, 2. In this case 7(a) is required up to 0(a 2 ), 
C(oi), £(ck) and .Fo(a) up to 0(a) together with the two-loop running of a in 7(a) and 
one-loop running of a in ((a). 

The functions entering the evolution equation can be determined by comparing eq. (jSJ) 
expanded in the coupling constant to the asymptotic, i.e. leading in M 2 /Q 2 , fixed order 
result for the form factor. To compute this fixed order asymptotic result we apply the 
expansion by regions approach formulated in [21] and discussed using characteristic two- 
loop examples in It consists of the following steps: (i) consider various regions of a loop 
momentum k and expand, in every region, the integrand in Taylor series with respect to the 
parameters that are there considered small, (ii) integrate the expanded integrand over the 
whole integration domain of the loop momenta, (iii) put to zero any scaleless integral. In 
step (ii) dimensional regularization j2H] with d = 4 — 2e space-time dimensions is used to 
handle the divergences. The following regions are relevant in the considered version (i) of 
the Sudakov limit |2Tj : 



hard (h) 

1- collinear (lc) 

2- collinear (2c) 

soft (s) 



k~Q, 

k+ ~ Q, k^ ~ M 2 /Q , k~M, 
k^ ~ Q, k + ~ M 2 /Q , k ~ M , 
k~M. (4) 



Here fc± = ko ± k$, k = (k\, £2). By k ~ Q, etc. we mean that any component of k^ is of 
order Q. In one loop this leads to the following decomposition [S] 

jrW = (A h + A c + A S )^ B , (5) 



At 1} = ^(-- + T (21n(g 2 )-3)-ln 2 (g 2 )+31n(g 2 ) + ^ 
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A« = C F ( ^ - i (2 ln(Q 2 ) - 4) + 2 ln(g 2 ) ln(M 2 ) - ln 2 (M 2 ) - 4 ln(M 2 ) - ^- + 4j , 

A« = C F (-1 + ln(M 2 ) + I) . (6) 

where = (A^ 2 — l)/(2iV) is the quadratic Casimir operator of the fundamental represen- 
tation of the SU(N) group and the subscript c denotes the contribution of both collinear 
regions. The 't Hooft scale [i has been dropped in the argument of the logarithms as well as 
the factor (47re _7B<E (/z 2 )) e per loop. For a perturbative function f(a) we define 
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The contribution of all the regions add up to the well known result 




On the other hand the one-loop form factor can be written as 

3=03 = (b m " l2 (w) + + (m ) ln (w) + F » a) ) F - ■ (9 > 

The expansion by regions is very efficient for determination of the parameters of the evolution 
equation. Indeed, it not necessary to compute the complete asymptotic result in order to 
obtain the functions parameterizing the logarithmic contributions. In the process of the 
scale separation through the expansion by regions the logarithmic contributions show up 
as the singularities of the contributions from different regions which cancel in the total 
result. Thus one can identify the regions relevant for determining a given parameter of the 
evolution equation and compute them separately up to the required accuracy. For example, 
the anomalous dimensions 7(a) and ((ce) are known to be independent on the infrared 
cutoff and are completely determined by the contribution from the hard loop momentum 
|18| H5] . If dimensional regularization is used for the infrared divergences of the hard loop 
momentum contribution, as in our approach, the anomalous dimensions 7(a) and ((ct) are 
given by the coefficients of the double and single poles of the hard contribution to the 
exponent (|HJ), respectively 01221 • On the contrary, the functions and F (a) fix the 
initial conditions for the evolution equation. They are not universal and depend on the 
infrared sector of the model. Furthermore, the values of and F (a) depend on the 
definition of the lower integration limits in eq. Q. To determine the function one has 
to know also the singularities of the collinear region contribution while i*o(a) requires the 
complete information on the contributions of all the regions. 

From the first line of eq. © we find the one-loop anomalous dimensions 

7 {1) = ~2C F , 

C (1) = 3CV- (10) 

With the above values of 7W and £W it is straightforward to obtain the one-loop result for 
the remaining functions 

= 0, (11) 

by comparing eqs. Q and (jSJ). Note that in the Born approximation Fq ^ = 1. 
A similar decomposition can be performed in two loops 

T {2) = (A hh + A hc + A cc + . . .) T B , (12) 
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Only the hard-hard part is now available [28]. However, this information is sufficient to 
determine the two-loop value Beside the running of the coupling constant, 7 < - 2 - ) is the 
only two-loop quantity we need for the NNLL approximation. It reads E 
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(2) 




20 



T F n f 



(13) 



for a defined in the MS scheme. Here Ca = N is the quadratic Casimir operator of the 
adjoint representation, Tp = 1/2 is the index of the fundamental representation and rif is 
the number of light (Dirac) fermions. 

Let us consider the two-loop corrections. The LL, NLL and NNLL approximations are 
given respectively by 



T (2) 
-r LL 



T {2) 
-> NLL 



T (2) 
■J I 



NNLL 



^(7 (2) + (C (1) -/5o)C (1) + Fo (1) 7 (1) )ln 2 (jP) 




(14) 
(15) 
(16) 



where /3o = 11CU/3 — 4Tpn//3 is the one-loop /3-function provided the normalization point 
of a is M. The two-loop running of the coupling constant in the leading order 7(a) starts 
to contribute in the three-loop NNLL approximation. 

The presence of a scalar particle in the fundamental representation with no Yukawa 
coupling to fermions leads only to a modification of the 7- and /3-f unctions. One scalar 
boson with the mass much less than Q gives the additional contribution of 16CVTf/9 to 7^ 2 - ) 
and the additional contribution of — Tp/3 to /3q- 

For the standard model inspired case of the SU{2) L and U{1) gauge groups with rif = 6 
and one charged scalar boson either in the fundamental representation of SU (2) or of the 
unit U(l) charge up to NNLL approximation we have 




(17) 



and 



-In 2 




\M 2 J 








V M 2 J 




+ 3 In 



52, 3 
— In , 

9 \M 2 
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respectively. The relatively small coefficient of the LL terms and the large coefficient of 
the NNLL terms in the form factor are clearly indicative of the importance of the NNLL 
corrections and, as we will see, reflect the general structure of the logarithmically enhanced 
electroweak corrections. 



3 The four fermion amplitude 

Let us now investigate the four-fermion scattering at fixed angles in the limit when all the 
invariant energy and momentum transfers of the process are far larger than the gauge boson 
mass, \s\ ~ \t\ ~ \u\ ^> M 2 . The analysis of the four fermion amplitude is complicated by 
the extra kinematical variable and the presence of different "color" and Lorentz structures. 
We adopt the following notation 

A x = ^°7/^i^°7/^3, 
All = ^2 L t a l^i L ^ L t a ^3 L , (19) 



etc. Here t a denotes the SU(N) generator, pi the momentum of the ith fermion and pi, p% 
are incoming, and p 2 , Pi outgoing momenta respectively. Hence t — (pi — p 4 ) 2 = — sa;_ and 
u = (p\ + ps) 2 = —sx + where x± — (1 ± cos6*)/2 and 9 is the angle between the spatial 
components of p\ and p^. The complete basis consists of four independent chiral amplitudes, 
each of them of two possible color structure. For the moment we consider a parity conserving 
theory, hence only two chiral amplitudes are not degenerate. The Born amplitude is given 
by 

A B = l -^-A x . (20) 
s 

The collinear divergences in the hard part of the virtual corrections and the corresponding 
"collinear" logarithms are known to factorize. They are responsible, in particular, for the 
double logarithmic contribution and depend only on the properties of the external on-shell 
particles but not on a specific process [13 HH H3 H3 UH 120] ■ This fact is especially clear if 
a physical (Coulomb or axial) gauge is used for the calculation. In this gauge the collinear 
divergences are present only in the self energy insertions to the external particles fUJ El 
ED] . Thus, for each fermion-antifermion pair of the four-fermion amplitude the collinear 
logarithms are the same as for the form factor T discussed in the previous section. Let us 
denote by A the amplitude with the collinear logarithms factored out. For convenience we 
separate from A all the corrections entering eq. © so that 

A=^(^) 2 A. (21) 



B 



The resulting amplitude A contains the logarithms of the "soft" nature corresponding to the 
soft divergences of the hard region contribution and the renormalization group logarithms. 
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It can be represented as a vector in the color/chiral basis and satisfies the following evolution 
equation [20112111231: 







; A = X (a(Q 2 ))A. 



(22) 
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where x( a ) is the matrix of the soft anomalous dimensions. Note that we do not include 
to eq. (J22)) the pure renormalization group logarithms which can be absorbed by fixing the 
normalization scale of g in the Born amplitude (J2(Jj) to be Q. The solution of eq. (J22|) reads 



Q 2, 
M 2 X 



(23) 



where Xi( a ) are eigenvalues of x( a ) an d Ah(c0 are Q-independent eigenvectors of x{ a ) which 
determine the initial conditions for the evolution equation at Q — M. Similar to the function 
F Q (a) they get contributions from all the regions while the matrix of the soft anomalous 
dimensions is given by the coefficients of the single pole of the hard region contribution to 
the exponent ([2"3~|) 1201 ■ Strictly speaking the matrices x( a (Q 2 )) f° r different values of 
Q do not commute and the solution is given by the path-ordered exponent [20] • The NLL 
approximation is given by the one-loop value of xi a ) while the NNLL approximation requires 
Aoi(a) up to 0(a) together with the one-loop running of a in x( a )- 

In one loop the elements of the matrix xi a ) do not depend on chirality and read [0] 



Xxx 



Xxl 



Xdx 
Xdl 



-2C A (In (x + ) + m) + 4 (c F - ^) In 



N 



x_ 



4 In 



0. 



x A 



X- 



(24) 



In higher orders the matrix xi a ) ma y De non- degenerate for the different chiral components 
of the basis. In the Abelian case, there are no different color amplitudes and there is only 
one anomalous dimension 



X 



(i) 



4 In 



x. 



(25) 



In terms of the functions introduced above the one-loop correction reads 



A® 



s 



7 v 




(26) 
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where A^ 



T-A {1) 



AS 1 '\q2 = m 2 has the following decomposition 



(i) 



A 



(1)A 



LL-^LL 



A 



A 



^0 LR-^-LR i * * * ■ 



(27) 



For the present two-loop analysis of the annihilation cross section only the real part of the 
coefficients A^ is needed, 



Re 
Re 
Re 
Re 



A (1)A 
A 
A 
A 



LL 

(1)A 
LR 

(l)d 
LL 

Li? 



7> 

TV 



Ca 
2 

CpT F 



85 



7T 



20 



N 
N 



N 

f\ x +i x ~ 



f(x-.,x+) + Cj 



85 2 

h 7T 2 

9 



20 



(28) 



f(x-,x + ) 



where 



f(x+, X- 



lnx_ 



X + i 2 

m x. 



X-\. 



xz 



(29) 



A scalar particle in the fundamental representation with no Yukawa coupling to fermions 
gives the additional contribution of — 87V/9 to the first two lines of eq. (J2EJ)- 

The two-loop correction is obtained by the direct generalization of the form factor anal- 
ysis. The only complication is related to the matrix structure of eq. 



A {2) 



LL 



A {2) 



NLL 



A {2) 



ig 2 (7 (1) ) ; 

s 2 

w 2 (Q 2 ) 

s 

ig 2 {Q 2 ) 



M 2 



A' 



NNLL 



In 4 

7 (2) + f2C (li - *nU' J ' 



Xxx 



2FiV ) + ^((4C (1) -/3o)xS+xir 

+x£!xSS)) ^ A + \ ((4C (1) - a.) x2 + xSSxS + x2x2) ^ + 7 (1) 4 X) 



(30) 
(31) 



xln 2 



M 2 



(32) 



The structure of the infrared singularities of the hard part of the two-loop corrections pre- 
sented here is in full agreement with the result of 30 which was confirmed by explicit 
calculation |3*T] . 

To illustrate the significance of the subleading contributions let us again discuss the 
standard model inspired example considered in the previous section. Having the result for 
the amplitudes it is straightforward to compute the one- and two-loop corrections to the 
total cross section of the four-fermion annihilation process using the standard formulae. For 
the annihilation process one has to make the analytical continuation of the above result to 
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the Minkowskian region of negative Q 2 = — s according to s + iO prescription. Although the 
above approximation is formally not valid for small angles 9 < Mj ^fs we can integrate the 
differential cross section over all angles to get a result with the logarithmic accuracy. In this 
way we obtain for the case of SU(2)l group 



and 



r(l) 



a® = 



80 

H T 

449 
IT 



In 

In 3 



-3 In' 



(m 2 ) 



26 



In 



(m 2 ) " 
(m 2 ) 



(m 2 ) 



25 



9 

/4855 



211 



+ 3vr' 



O'B , 



9 



+ 



3vr' 



37tt 2 



In' 



(m 2 ) 



Ob , 



2 



799 37tt 2 



vm 2 ; 



6 



\M 2 J 



9 



cb , 



In' 



(m 2 ) 



0"B , 



(33) 



(34) 



for the initial and final state fermions of the same or opposite isospin, respectively. Here Ob 
is the Born cross section with the MS couplings constant normalized at the scale y/s. 
For the U(l) group we have 



cr ' = 



-4 In' 



(m 2 ) 



+ 12 In 



r ( 2 ) 



532 
IT 



lu 



(m 2 ) ' 

VM 2 / 



382 4vr 2 



9 



+ 



1142 167r 2 '| h| 



2 f — ) 

VM 2 / 



(35) 



Similar to the form factor, we observe a relatively small coefficient of the LL terms and a 
large coefficient of the NNLL terms. 



4 NNL logarithms in electroweak processes at high en- 
ergies 

We are interested in the process /'/' f f ■ In the Born approximation, its amplitude is of 
the following form 

frw +t 2 «> l£ p-)M'j, pel 

6 I,J=L,R V ^ 7 

where 

M'l = fMJflpfj , (37) 

tvK = tan6*vK with being the Weinberg angle and Tf (Yf) is the isospin (hypercharge) of 
the fermion which depends on the fermion chirality. 

To analyze the electroweak correction to the above process we use the approximation 
with the W and Z bosons of the same mass M, the Higgs boson of the mass Mh ~ M and 
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massless quarks and leptons. A fictitious photon mass A has to be introduced to regularize 
the infrared divergences. We insert the mass into the gauge boson propagators "by hand" to 
investigate the leading in s _1 behavior of the amplitudes, leaving aside the Higgs mechanism 
of the gauge boson mass generation. This approach is gauge invariant as far as power 
unsuppressed terms are considered. The NNLL approximation is not sensitive to the fine 
details of the mass generation because we need only the hard part of the potentially dangerous 
self-energy insertion to the gauge boson propagator. Indeed, the only effect of the virtual 
Higgs boson is the modification of the functions (3q, 7^ and A^ x . The function as well 
as the running of the coupling constant in 7' 1 ' and are determined by the singularities 
of the hard part of the corrections. On the other hand, the vacuum polarization of the 
off-shell gauge boson in the Born amplitude contributing to A^ x is infrared safe and can be 
computed in the massless approximation, i.e. in the leading order in s _1 it receives only the 
contribution of the hard region. The only effect of the Higgs mechanism is that we have two 
different masses M z and M w = cos 9wM Z - Since cos 6w ~ 1 we neglect this difference in our 
calculation. The correction due to the heavy gauge boson mass splitting will be discussed at 
the end of the section. 

Let us first consider the equal mass case A = M, where we can work in terms of the fields 
of unbroken phase. In the massless quark approximation the Higgs boson couples only to 
the gauge field. Therefore the result of Sects. 2 and 3 for the SU(2)l gauge group with the 
coupling g and the U(l) gauge group with the coupling twQ can be directly applied to the 
electroweak processes. For the standard electroweak model with one charged Higgs doublet 
one has to replace in the above expressions rif — > 2N g + 1/4 for SU(2)l /^-function, Tprif — > 
5N g /3 + 1/8 for U(l) ^-function, n f -> 2N g + 2/5 for SU(2) L 7 (2) and coefficients, and 
TpUf — > 5N g /3 + 1/5 for U(l) 7*- 2 - 1 and Aq coefficients, with N g = 3 being the number of 
generations. For example we have 

7 (2) 

A) 

for SU(2) L and 

7 (2) 

Po 

for 17(1). 

The result for the amplitudes is obtained by projecting on a relevant initial/final state 
with the proper assignment of isospin/hypercharge. For example, the projection of the ba- 
sis ()19|) on the states corresponding to the neutral current processes reads A)j -> TfTftAfJ 
Ajj — > AjJ . The only complication in combinatorics is related to the fact that now we are 



10 Ar 65 2 

— N„ + 7T 2 

3 9 3 
4 Ar 43 

-3 N - + T' 



200 Ar 8 
27 9 9 

20 ,r 1 

■— N n - - 



V 2 

l W A 



(38) 



(39) 
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having different gauge groups for the fermions of different chirality. In particular, the double 
logarithmic approximation is given by the exponential factor 



exp 



where 



-\T f (T f + r 



HQ 2 



Y 2 



HQ 2 



r 



16vr 2 



In 2 




(40) 



(41) 



and 27(27 + 1) = C F . 

The photon is, however, massless and the corresponding infrared divergent contributions 
should be accompanied by the real soft photon radiation integrated to some resolution energy 
u res to get an infrared safe cross section independent on an auxiliary photon mass. In 
practice, the resolution energy is much less than the W (Z) boson mass and the massive gauge 
bosons are supposed to be detected as separate particles. To study the virtual corrections in 
the limit of the vanishing photon mass we follow a general approach of the infrared evolution 
equations developed in jB] (see also references therein). It is convenient to use the auxiliary 
photon mass A as a variable of the infrared evolution equation below the electroweak scale 
M. The dependence of the virtual corrections on A in the limit A <C M is canceled by the 
contribution of the real soft photon radiation. For u res <C M, the soft photon emission is of 
the pure QED nature. Therefore, the kernel of the infrared evolution equation which governs 
the A dependence of the virtual corrections to the amplitudes is Abelian. This dependence is 
given by the QED factor U which can be directly obtained from the general formulae given 
above: 



U 



U (a e )exp 




a e (A 2 
Atx 



Q) + <tr "(I (<JJ + Q 2 r) + 



- (3 {Q) + Q)) + iQ,Q r In gi)) In (J) + £ {Q) + Q)) ^N s 
+ 0{<4)\, (42) 



where a e is the MS QED coupling constant and we use the following expressions for the 
QED functions 



(i) 



c 



3Q 



f ■ 



4Qf>Qf In I 

N a , 

9 

27 9 f 



X j 



X. 



(43) 
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The expressions for /3q and % 2 > can be obtained by substituting TpUf — > SN g /3 to the general 
formulae. The coefficient Uo(a e ) in eq. ()42|) is a two-component vector in the chiral basis. 

In a full analogy with the renormalization group all the information on the non-Abelian 
gauge dynamics above the electroweak scale up to power suppressed contributions is con- 
tained in the initial condition for this Abelian infrared evolution equation at the point A = M. 
To fix a relevant initial condition for the evolution in A below the electroweak scale one has 
to subtract the QED virtual correction (}4~2*|) computed with the photon of the mass M from 
the complete result with A = M jH]. This leads, in particular, to the modification of the 
function 7(a) so that the double logarithmic exponential factor becomes 



exp 



~[T,(T f + l) 



4j - W + (/ ~ /')) HQ 2 



(44) 



where sw = sm8\y In the NNLL approximation after the subtraction we get 



7 



(2) 



20, T 130 2tt 2 

— N a H 

9 9 9 3 



) T d T, + 1) - + 



160 



27 



Y 2 
f 4 -L 



(45) 



A similar subtraction should be done for the parameters Q 1 ' and which take the form 



and 



(46) 



Xxl 



Xdd 



-4 (In + ztt) + (4l7,y> - As^Qf'Qf + 2) In 



4 In 





t^YfYf - AswQfQf) In I -pi . 



(47) 



x — j 



For I or J = R the matrix is reduced to 



(i) 



41>1) - te^QpQ, 



In 



X - 



X- 



(48) 



At the same time we have some freedom in the definition of the coefficients Fo(a), Ao(a) 
and U (a e ). If we use the one-loop normalization condition U\q2 =M 2 = 1, then Uq 



(i) 







13 



and no QED subtraction is necessary for F (1) and A^. In this case the QED factor U is 
universal and has no matrix structure. To summarize, we have two evolution equations and 
corresponding initial conditions. The coefficients Fq{q) and Aq{q) give the initial condition 
for the hard evolution of the amplitudes in Q at Q = M while the above subtraction of the 
QED contribution gives the initial condition for the infrared evolution of the amplitudes in 
A at A = M. 

The result for the n-loop correction to the amplitude (J36|) can be decomposed as 

A™ = a£1 + aP ll + aW nll + .... (49) 

Explicit expressions for A LL and A NLL can be found, for example, in In the NNLL 
approximation one has to take into account also the one-loop constant contribution corre- 
sponding to Fq 1 ^ and Ajp terms of eq. (j2Tlj) . It reads 



a rt-NNLL 



(50) 



* I,J=L,R 



+ ( 2T}Tf + t 2 w 



Y f' Y f 



2vr 2 



7)07) 



i) + 4^f + (/«-/') 



T*Tf + 4 



2 Y f' Y f 

w 



+ 



4 y w 4 

f(x-,x + ) (5 IR 6j L + S IL 6 JR ) 

2tt 2 



f(x+, x-) (5 IR 5 JR + 5 IL 5 JL ) 

2(1 4\ „ „ /100 AT 4 
9 9 9 J f f V 27 9 9 



,4 *>*7 



where a = (7 2 /167r 2 and we keep only the real part of Aq . 

Let us consider the two-loop corrections. The two-loop LL corrections to the chiral 
amplitudes were obtained in jB] 



s 



I,J=L,R 



(51) 



~ (T f (T f + 1) + t\ Y l - s 2 w Q) + (/ <- /')' 



TpTf + t 2 w - 



Y f Y f 



4 



and the two-loop NLL corrections with the exception of the trivial corrections proportional 
to (3q can be found in P 

,2 



a 2 A 



(2) 
NLL 



* I,J=L,R 



T f (T f + 1) + 4-2 _ ^Qj + (/<_> /') 



(52) 



x 1 3 
+ - 
+ ln 



T /( T / 
4 (In (a^ 



i) + 4? -swQ) + (/ W) 



rp3 rp3 . ,2 f/j/ 



+ «r) + In (2 + t 2 w Y f Y f )j Tf,Tj + | In ^ 



T"i3 rp3 | ,2 f/j/ 



L(Q 2 )/(Q 2 )A{/ 
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where 



KQ 2 ) 



167T 2 11 U 2 



(53) 



with 5n = 1 for I = L and zero otherwise. The second line of eq. (J52j) corresponds to the 
term of eq. (J3~T]) while the third and forth lines correspond to the terms in eq. ([31)1 . 
A part of the (3q NLL corrections is absorbed by choosing the normalization point of the 
coupling constants in eq. (J31j|) to be Q. The rest is due to the running of the coupling 
constant in the double logarithmic integral and corresponds to the /3q term in eq. (|3ip. It is 
of the form 



n 2 4 < 2 ) I 
" ^nllIPo 



1 

3 

x 



s 
43 



E 

I,J=L,R 



3 9 6 



Tf(T f + l) + 



9 9 



KQ 2 )L(Q 2 )M'/ 



(54) 



32 



4 



9 



i j i 



provided the normalization point of the coupling constants is M with the exception of the 
coupling constants entering the Born amplitude (j3T?j) normalized at the scale Q. 
Let us consider the NNLL contribution. For convenience we split it in four parts 



.4 



(2) 

NNLL 



A /t( 2 ) i A A {Z> \ \'~' _i_ A A^> 

I NNLL + 2 NNLL ~+~ ^ A NNLL ~+~ ^ A NNLL 



| (2) 



(2) 



| (2) 



(55) 



where the trivial 0k renormalization group logarithms which can be absorbed into the running 



of the coupling constants in eq. (fBTIj) are not included. The correction A 1 A) v j Vii correspond- 
ing to the 7W, C (1)2 , /?oC (1) and F (1) 7 (1) terms of eq. © is 



^ 2 A 4 (2) 



W_ 

s 



E 

I,J=L,R 



(56) 



20 Ar 130 
9 9 9 



2tt 2 



Tf(T f + r 



/100 Af 4 



27 



+(/ - /') 



+- 



41 



-Nn - 

3 9 
2tt 2 



9 

+ 2 
43 



T/Cr, + i) + 4-2 - 4qJ + (/ <-> /') 



r 7 (T / + i) 



4 

20 Ar 1 

9 9 6 



y 



32 



2 3 
+ (/ - /' 



Y 2 



T f (Tf + l)+t 2 w ^--s 2 w Q} + (f~f) 



Tf(T f + l)+t 2 w ^l 



TfiTf + 1 



Y r Yf 



w 



i\Q 2 )Ag. 
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The correction A 2 v4^]y LL corresponding to the C^X^ an d A)X^ terms of eq. (jS2} reads 

.2 



a 2 A 2 A 



(2) 

NNLL 



E 



(57) 



I,J=L,R 



-4 (In (x_ 



in) + In 



X 



2 + t z w Y f ,Y f )\nT* + ^\n 



- + 3 (Tf(T f + 1) + - s^Q 2 . + (/ ~ /') 



+4 In — 



^rAT + — U 



9 



12 



y2 

4 



x + 



X. 



5tt,5 



il<jjl 



Y N 9 s wQf>Qf)+3[t 



a Yf'Yf 

■W A 



swQrQf 



x (r f {T f + 1) + - s 2 w Q) + (/<-> /') 



The correction A 3 A^ NLL corresponding to the terms of eq. (}32|) reads 



rp3 rp3 i ,2 £/j7 



l 2 (Q 2 )4J 



^■i^NNLL — E 

S I,J=L,R 



(58) 



lH-4(ln( ; r + )+i7r)+l n ^( 

-4 (In (x + ) + ztt) + 2 In ^ (l + t 2 w Y f ,Y f - 4^^) 



'2 + tj^fyy, - 4*iyQ r Q,) I T f 3 ,T f 3 + In 



X 



\5 IL 5j L + e w Y s ,Y f T}Tj 



+ In 



2 I X + 



x_ 



:\[nT* + t 2 w ^-t5 IL 5 JL 



+ {t 2 w Y r Y f -As 2 w Q r Q f ) t 2 w - 



2 2 *>!7 



The correction A4 J 4^] Vii corresponding to the real part of the j^A^ term of eq. (|32|) reads 
a 2 A 4 ^gU = -^ E (59) 



I,J=L,R 



3^3 , + 2 J '/' J '/\,2 Yf'Yf 



(r f ( Tf + 1) + 4^ _ ^ + (/ /0 j I ( 27?i7? + 4 ^ 



X 



ac_) (6 IR 5 JR + S IL 5 JL ) - /(ar_, x+) (5 /i? 5 JL + 5 IL S JR ) 



20 4\ o , 
Y 9 + 9) T /' T / 



100 A7 4, 
27 9 9 



4 *W 
VK a 



1. , 170 n 2 



7/,T| + x_)5 /L 5 JL } / 2 (g 2 )A{7 . 



With the expression for the chiral amplitudes at hand, we can compute the leading and 
subleading logarithmic corrections to the basic observables for e + e~ — > ff. 
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In the NNLL approximation one has to take into account also the effect of analytical 
continuation to the physical positive real value of the invariant s. For the annihilation 
processes it is more natural to normalize the QED factor at the Minkowskian point s = M 2 
to U\ s= m2 = 1 so that after the expansion in a e it reads 



U = {l-^ (Q} + Qf) ^ (£) - ((3 + 2^) (Q) + Q)) + 4 Q f Q f , In ( 



x In ( -^r 



A 



+ 0{al) ) , (60) 



Let us consider the total cross sections of the quark- antiquark//i + /i~ production in the 
e + e~ annihilation. The LL, NLL and NNLL corrections to the cross sections to one and two 
loops read 

Rqq = l-1.66L(s)+ 5.3H(s)- 8.36a + 1.93 L 2 (s) - 10.59 L(s)/(s) + 31.40/ 2 (s), 
R q g = 1-2.18 L(s) + 20.58 l(s) - 34.02 a + 2.79 L 2 (s) - 51.04 L(s)l(s) + 309.34 l 2 (s) , 
Rp +Ii - = 1 - 1.39 L(s) + 10.12 l(s) - 20.61 a + 1.42 L 2 (s) - 19.81 L(s)l(s) + 107.03 l 2 (s) , 

(61) 



where Q = u,c,t, q = d,s,b, Rqq = a/<jB(e + e~ — > QQ) and so on. The MS couplings in 
the Born cross section are normalized at y/s. Numerically, we have L(s) = 0.07 (0.11) and 
l(s) = 0.014 (0.017) for a/s = 1 TeV and 2 TeV respectively. Here M = Mw has been chosen 
for the infrared cutoff and a = 2.69 ■ 10~ 3 , = 0.231 for the MS couplings normalized 
at the gauge boson mass. The small difference between the two-loop NLL coefficients in 
eq. (JoTj) and the result of P is due to the /?o contribution (p>Ij) . 

To get the infrared safe result for the semi-inclusive cross sections one has to add to the 
expressions given above the standard QED corrections due to the soft photon emission and 
the pure QED virtual correction which is determined for massless or light fermions of the 
mass mj < A < M by eqs. (}4*2*)) . ()60|). To derive the QED factor for A far less than the 
fermion mass A rrif M one has to change the kernel of the infrared evolution equation 
and match the new solution to eq. ()42|) at the point A = rrif. The sum of the real and 
virtual QED corrections depends on s, u res and on the initial/final fermion masses but not 
on Mzyv- Note that our analysis implies the resolution energy for the real photon emission 
to be smaller than the heavy boson mass. If the resolution energy exceeds Mz,w the analysis 
is more complicated due to the fact that the radiation of real photons is not of Poisson type 
because of its non-Abelian SU(2)l component 0. In the case of the quark-antiquark final 
state the strong interaction also produces the logarithmically growing terms. They can be 
read off the results of Section 1 for the form factor. For massless quarks the complete 0(a 2 ) 
corrections including the bremsstrahlung effects can be found in |2*5] . 

For completeness we give a numerical estimate of corrections to the cross section asym- 
metries. In the case of the forward-backward asymmetry A FB (the difference of the cross 
section averaged over forward and backward semispheres with respect to the electron beam 
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direction divided by the total cross section) we get 

1 - 0.09L(s) - 1.23Z(s) + 1.47a + 0.12 L 2 (s) + 0.64 L(s)l(s) - lA0l 2 (s), 
1 - O.UL(s) + 7.15 l(s) - 10.43 a + 0.02 L 2 (s) - 1.31 L{s)l{s) - 33.46 l 2 (s) , 
1 - 0.04 L(s) + 5.49 /(s) - 14.03 a + 0.27 L 2 (s) - 6.32 L{s)l{s) + 21.01 l 2 (s) , 

(62) 

where R FB = A FB /Ag B . For the left-right asymmetry A LR (the difference of the cross 
sections of the left and right particles production divided by the total cross section) we 
obtain in the same notation 

2.34L(s)+ 8.98 l(s) - 5.73 a - 0.46 L 2 (s) + 7 A3 L(s)l(s) - 18.59 l 2 (s) , 
1.12 L(s) + 11.86 l{s) - 15.83 a - 0.81 L 2 (s) + 17.74 L(s)l(s) - 127.05 Z 2 (s) , 
13.24 L(s) + 113.77 l(s) - 139.94 a - 0.79 L 2 (s) + 23.34 L(s)Z(s) - 155.36 l 2 (s) . 

(63) 

Finally, for the left-right asymmetry A LR (the difference of the cross sections for the left and 
right initial state particles divided by the total cross section) which differs from A LR for the 
quark- ant iquark final state we have 

1- 2.75L(s)+ 10.07 /(s)- 9.02 a - 0.91 L 2 (s) + 10.80 L(s)l(s) - 32.10 l 2 (s) , 
1 - 1.07 L(s) + 11.56 /(s) - 15.60 a - 0.77 L 2 (s) + 16.78 L(s)l(s) - 121.56 l 2 (s) . 

(64) 

In the 1 — 2 TeV region the two-loop LL, NLL and NNLL corrections to the cross sections 
can be as large as 1 — 4%, 5 — 10%, and 5 — 9% respectively. However, we observe a significant 
cancellation between different terms and the sum of the known two-loop corrections amounts 
of approximately 1 — 2%. The sum of the two-loop correction to the asymmetries is even 
smaller and does not exceed 1% level with the exception of the • For this quantity the 

relatively large corrections are the consequence of the numerically small Born approximation. 

Let us discuss the accuracy of our result. At TeV energies the LL, NLL and NNLL 
corrections of eqs. provide asymptotic expressions for the cross section and the 

asymmetries to one and two loops. The complete one-loop corrections are known exactly 
(see (221 f° r t ne most general result) and we have included the dominant one- loop terms in 
eqs. (|61 )) -([64 )) to demonstrate the structure of the expansion rather than for precise numerical 
estimates. For physical applications, the mass difference between the W and the Z gauge 
boson, power suppressed terms and also top quark mass effects can be important. The effect 
of the W and Z gauge boson mass difference on the coefficients of the NLL and NNLL terms 
is suppressed as [Mz ~ My/)/M ~ 0.1 while the leading power corrections can be as large 
as M 2 /s < 0.01. Thus, except for the production of third generation quarks, the above 
expressions approximate the exact one-loop result with 1% accuracy in the TeV region. At 
the same time both effects can be neglected in two-loop approximation. Therefore, the only 
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essential deviation of the complete two-loop NLL and NNLL result from eqs. (|61| ) -([64 p for 
the production of the third generation quarks is due to the large top quark Yukawa coupling. 
Numerically, the corresponding corrections can be as important as the generic non- Yukawa 
ones. 

Finally, let us emphasize that the angular dependent NLL and NNLL terms are quite 
important for the cross section and dominate in particular the forward-backward asymmetry. 

5 Summary 

In the present paper we employed the evolution equation approach to analyze the high energy 
asymptotic behavior of the four-fermion amplitudes in the non-Abelian gauge models. The 
results were used to compute the NNLL electroweak corrections to the neutral current four- 
fermion processes at high energy in the massless quark approximation to all orders in the 
coupling constants. We have shown the NNLL approximation to be insensitive to the details 
of the gauge boson mass generation as well as to the Higgs boson mass and self-coupling. 

We have calculated the explicit expressions for the one- and two-loop terms which saturate 
the NNLL corrections to the basic observables in the TeV region. In general, the two- 
loop NLL and NNLL corrections exceed the LL contribution in the TeV region due to the 
numerically small coefficient in front of the double logarithmic terms. Hence the truly 
asymptotic behavior sets in only at an significantly higher energy. At the same time the 
two-loop NNLL corrections are numerically of the same magnitude but slightly smaller than 
the NLL ones, both being in the range of 1%-10%. This could be considered as a signal of 
convergence of the logarithmic expansion at TeV energies. Indeed, the two-loop coefficients 
in front of In 2 (s/M 2 ) is (a few units) x a 2 . This is not an unusually large value in a non- 
Euclidean regime where the expansion parameter is a, rather then a/(47i~) as can be seen in 
eqs. ([3~3*)) - ([3l3)) (see also [221 I2H EH] ) • Moreover, we have observed a significant cancellation 
between the two-loop LL, NLL and NNLL terms. As a result of this cancellation the sum of 
these two-loop corrections to the cross sections is of order 1 — 2% for all the processes. 

Thus, if we assume no further growth of the coefficient for the single logarithmic and 
non-logarithmic two-loop terms and the observed pattern of cancellation to hold, we would 
argue that our NNLL result approximates the exact cross sections with 1% accuracy. The 
accuracy is less for the production of third generation quarks where we cannot neglect the 
large Yukawa coupling to the Higgs boson that modifies the NLL and NNLL terms in our 
formulae. 
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